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Abstract. Based on correction from instanton-gluon interference to correlation function, the properties of 
the 0 ^ pseudoscalar glueball is investigated in a family of finite-width Gaussian sum rules. In the frame¬ 
work of semiclassical expansion for quantum chromodynamics (QCD) in the instanton liquid background, 
the contribution arising from the the interference between instantons and the quantum gluon fields is cal¬ 
culated, and included in the correlation function together with pure-classical contribution from instantons 
and the perturbative one. The interference contribution is turned to be gauge-invariant, free of infrared 
divergence, and has a great role to restore the positivity of the spectra of the full correlation function. 

The negligible contribution from vacuum condensates is excluded in our correlation function to avoid the 
double counting. Instead of the usual zero-width approximation for the resonances, the usual Breit-Wigner 
form with a suitable threshold behavior for the spectral function of the finite-width resonances is adopted. 

A consistency between the subtracted and unsubtracted sum rules is very well justified. The values of the 
mass, decay width and coupling constants for the 0 ^ resonance in which the glueball fraction is dominant 
are obtained, and agree with the phenomenological analysis. 
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1 Introduction 

A significant issue in quantum chromodynamics (QCD) 
is to seek for the signal of the existence of glueballs. Be¬ 
cause glueballs are bound sates composed of only gluons 
in the quarkless world, such signal may give a unique in¬ 
sight into the non-Abelian dynamics of QCD. Theoretical 
investigations including lattice simulations model 

researches [HEIE] and sum rule analyses [TlIRlIlIl ITTIlimfT^ 
have been going on for a long time, but no decisive evi¬ 
dence of the existence of glueballs has been confirmed by 
experimental research up to now [nin]. Further investi¬ 
gation on glueballs still makes sense. 

One of the obstacles in theoretical researches of glue¬ 
balls is that non-perturbative dynamics of QCD, which 
is responsible for the formation of hadrons, is difficult 
to handle, and the QCD vacuum is recognized to be a 
medium with complicated structure, and may impact greatly 
on the attributes of hadrons. In particular, the tunneling 
effect between the degenerate vacua of QCD should be 
taken into account. In the leading order, this effect is de¬ 
scribed by instantons and shown to be of great 

significance in generating the properties of the unusual 
hadrons, glueballs. Moreover, the glueball may be mixed 
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with usual mesons of the same quantum numbers, making 
the identification of the glueball more complicated UMH- 

Instantons, as the strong topological fluctuations of 
gluon fields in QCD, are widely believed to play an impor¬ 
tant role in the physics of the strong interaction (for re¬ 
views see mm)- In particular, instantons provide mech¬ 
anisms for the violation of both U{1)a and chiral sym¬ 
metry in QCD, and may therefore be important in deter¬ 
mining hadron masses and in the resolution of the famous 
17(1)A problem. Furthermore, it was recently shown that 
instantons persist through the deconfinement transition, 
so that instanton-induced interactions between quarks and 
gluons may underlie the unusual properties of the so called 
strongly coupled quark-gluon plasma recently discovered 
at RHIC m\- 

In the instanton liquid model, a narrow sense describes 
the QCD vacuum as a sum of independent instantons 
with radius p = (600MeV)“^ and effective density n = 
(200MeV)"'^ [10]. This model avoid the infrared problem 
caused by an infinite instanton density in the diluted gas 
model. The correctness of the instanton liquid model is 
still being intensively investigated. So far the model is 
essentially justified by its phenomenological success. The 
most important predictions are probably the breaking of 
the chiral symmetry (SBCS) in the axial triplet channel 
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mm and the absence of Goldstone bosons in the axial 
singlet channel. 

The non-perturbative effects of QCD is commonly casted 
into the description of non-perturbative vacuum, such as 
quark and gluon condensates and the instanton configu¬ 
rations. When contributions from both condensates and 
the instanton are included into the correlation function 
together, it leads to a double counting in a sense that 
they can be two alternative ways for parametrizing the 
nonperturbative vacuum mm- Moreover, the instanton 
distribution is closely connected with the vacuum conden¬ 
sates since the mean size and density of instantons can be 
deduced from the quark and gluon condensates, and con¬ 
versely, the values of condensates can be reproduced from 
instanton distribution . The contributions of 

instanton and those of the condensates should be equiv¬ 
alent from each other in some case, and thus including 
both contributions at the same time will cause the dou¬ 
ble counting problem. To avoid it, many authors invoke 
some special techniques in dealing with the two contribu¬ 
tions This issue is, however, really not set¬ 

tled. Fortunately, the contribution of condensates to the 
correlation function in the glueball channels is unusually 
weak as demonstrated in this work and by many other 
authors [33l[34]. Therefore, it is assumed that the conden¬ 
sate contributions can be understood as a small fraction 
of the corresponding instanton one in the local limit. To 
prevent thoroughly from the problem of double counting, 
we choose to work in the instanton vacuum model of QCD, 
and carry out a semi-classical expansion in instanton back¬ 
ground fields as suggested in our previous works to ana¬ 
lyze the properties of the lowest 0+“*' scalar glueball [53 
156] and 0 pseudoscalar one [571158] . where the correla¬ 
tion function of the glueball currents are calculated by just 
including the contributions from the pure instantons, the 
pure quantum gluons, and the interference between both, 
instead of working with both instantons and condensates 
at the same time. 

Zhang and Steele [55] have reported a disparity be¬ 
tween their Laplace sum rule and Gaussian sum rule for 
the pseudoscalar glueball. In this reference, the optimized 
parameters of pseudoscalar glueball have been obtained 
from the Gaussian sum rule, while it fails to achieve a sat¬ 
isfying Laplace sum rule. This result seems strange, since 
both Laplace and Gaussian sum rules are derived from the 
same underlying dynamical theory, and should, at least, 
be approximately consistent. It reflects the inconsistency 
for including both condensate and instanton contributions 
in the correlation function and disregarding the important 
contribution from the interaction between instantons and 
their quantum counterparts at the same time. 

An other serious problem in the 0 ^ glueball sum 
rule approach is that the fundamental spectral positivity 
bound is violated when including the strong repulsive pure 
instanton contribution, and as a consequence, the signal 
for the pseudoscalar glueball disappears [55]. To cure this 
pathology of positivity violation, the topological charge 
screening effect in the QGD vacuum is added to the cor¬ 
relation function, and a suitable instanton size distribu¬ 


tion is taken into account [4()l [7]. However, as comparing 
with the interference contribution, which we have recalcu¬ 
lated in this paper, and the pure perturbative one in the 
considered energy region, the topological screening effect 
turns out to be negligible; and the pathology of positiv¬ 
ity violation disappears when including the interference 
contribution (s. below). 

Phenomenologically, the identification of the pseudoscalar 
glueball has been a matter of debate since the Mark II ex¬ 
periment proposed glueball candidates [H]. Later, in the 
mass region of the first radial excitation of the rj and rj' 
mesons, a supernumerous candidate, the ?7(1405) has been 
observed. It turns out to be clear that 77(1405) is allowed 
as glueball dominated state mixed with isoscalar qq states 
due to its behavior in production and decays, namely, it 
has comparably large branching ratios in the J/ip radiative 
decay, but not been observed in 77 collisions [421IT51I45] . A 
review on the experimental status of the 7(1405) is given in 
Ref. [T5]. However, this state lies considerably lower than 
the theoretical expectations: the lattice QCD predictions 
suggest a glueball around 2.5GeV @3113; the mass scale 
of the pseudoscalar glueball obtained in the QCD sum rule 

approach is above 2GeV [5711551115113 ■ On the other hand, 

there are attractive arguments for the approximately de¬ 
generate in mass for the scalar and pseudoscalar glueballs 
@5] - and even the scenario that a pseudoscalar glueball 
may be lower in mass than the scalar one is recently dis¬ 
cussed in Ref. |47j . The possibly non-vanishing gluonium 
content of the ground state 7 and 7 ' mesons is discussed 
in @51155115511X5] . Up to now, only a topological model of 
the glueball as closed flux tube @3 predicts degeneracy of 
the O'*''*' and 0 glueball masses and admits the region 
1.3-1.5GeV. 

To face a more realistic phenomenological situation, 
we now reexamine the correlation function in the instan¬ 
ton liquid vacuum, and include all the resonances below 
and near the 7(1405) into the finite-width spectral func¬ 
tion, and then achieve a series of results in traditional 
Gaussian sum rule analyses which are consistent with the 
phenomenology. On the other hand, as a crosscheck, these 
results are almost same with the Laplace ones (paper is 
submitted), because both Laplace and Gaussian sum rules 
are derived from the same underlying dynamical theory. 
The paper is organized as follows: In the second section, 
we present systematically the calculation of the contri¬ 
bution to the correlation function due to the interference 
between instantons and quantum gluons. It is this inter¬ 
ference between instantons and quantum gluons to serve 
to be a mechanism to keep the positivity of the spectral 
function in contrast with the so-called topological charge 
screening effect stressed in Ref. [7] . The effect of the contri¬ 
bution of topological charge screening is found to be negli¬ 
gible as comparing with the interference one in the consid¬ 
ered energy region. The spectral function is constructed 
in the similar way as in case of the 0'''+ glueball [551I5I] 
with a suitable threshold behavior in the third section. 

In the fourth section, a family of Gaussian sum rules are 
constructed. The numerical simulations are carried out in 
the fifth section, and the results are consistent with vari- 
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ous Gaussian sum rules and in accordance with the phe¬ 
nomenology. Finally, the main conclusions are given, and 
a discussion of some interesting issues is open. 


2 Correlation function 


We are working in Euclidean QCD. The pseudoscalar 
glueball current is defined as 

Opix) = (1) 

where as is the strong coupling constant, G'^^^{x) is the 
gluon field strength tensor with the color index a and 
Lorentz indices fi and and 

Gfj,u,a{x) = l^^livpcrG p^ix) (2) 

is the dual of The current Op{x) is a Lorentz- 

irreducible, gauge-invariant and local composite operator 
with the lowest dimension, and renormalization group in¬ 
variant at least to the leading order as in the quark¬ 
less world. It is noticed that the current Op{x) is anti- 
hermitian due to the involving of the imaginary time, 
while its analytic continuation to the Minkowskian space- 
time is hermitian. The QCD correlation function is defined 
as 

n{q^) = J d'^xe^’^-^n\T{Opix)olio)}\n). (3) 


where Oj, is the hermitian conjugation of Op. The advan¬ 
tages to use the hermitian conjugation in definition are: 
first, the spectral functions both in Euclidean space-time 
and in its analytic continuation into Minkowskian space- 
time are, in principle, positively definite; and second, the 
relationship between the correlation functions both in Eu¬ 
clidean and Minkowskian formulations becomes very sim¬ 
ple as 


nE{Q^ = q^)^nM{Q^ = -q^), (4) 

because the overall minus sign arising from the analytic 
continuation due to {epupaip.i^pa)E = is 

just canceled with another minus sign arising from the 
mentioned different hermiticity of the pseudoscalar glue- 
ball current. Therefore, the expressions for JJEiQ'^) and 
are, in fact, the same function of Q^. 

In the framework of semiclassical expansion, the glue 
potential field B{x) can be decomposed into a summation 
of the classical instanton A and the corresponding quan¬ 
tum gluon field a as 


(5) 


Consequently, the pure-glue Euclidean action can be ex¬ 
pressed as 

5[i3] =So-l |l[A + a] + 

= So-lJd^x {a“ + 2gr>>^Fl 

‘^9 fabc^fib^ 


_ (1 _ ^ 


^9 fabcO’fibG' yc fadeO'fid^ V 


( 6 ) 


where Sq = %Td j is the one-instanton contribution to 
the action, Fp^a is the instanton field strength tensor 

Bpu.a^Ad^ — dpAif fi difAp si F gs fabcAp^i/Au Q^ (7) 

and Dp^{A) the covariant derivative associated with the 
classical instanton field 


D;\A) = dpSab + gfacbA'p 
In addition, the background field gauge 
Df^{A)ap = 0 


( 8 ) 

(9) 


is used with ^ being the corresponding gauge parameter, 
and certainly, the corresponding Faddeev-Popov ghosts 
according to the standard rule should be added to restore 
the unitarity. We note here that the structure constants 
fabc should be understood as eabc when any one of the 
color-indices a,b and c is associated with an instanton field 
due to the property of the closure of any group. 

According to the decomposition ([5|), the correlation 
function 77 splits into three parts, namely the pure classi¬ 
cal part, the pure quantum part and the interference part 
in the leading order 

77QCD(g2^ = 77('=9 (q 2)+ 7j(qu)(Q2)+i7M)(Q2^_ 


where the superscript indicates that it is calculated in the 
underlying dynamical theory, QCD. It is important to note 
that every part in rhs of (HUD is gauge-invariant because 
the decomposition ®, in principle, has no impact on the 
gauge-invariance of the correlation function. The pure in- 
staton contribution 77('^9 (^Q^) and the perturbative contri¬ 
bution 7T(‘1")(Q2) 

up to three-loop level in the chiral limit 
of QCD are shown in Eqs. o and dSU respectively in 
Appendix [A] The contribution of the topological charge 
screening which is not included in rhs of m is given in 
Appendix |B] 

One of our main tasks in this work is to calculate the 
contribution 7T^'"*)((5^) in (0, which is arising from the 
interference between the classical instantons and quantum 
gluons in the framework of the semi-classical expansion 
for QCD with the instanton background. After imposing 
the background covariant Feynman gauge (C = 1) for the 
quantum gluon fields, we are still free to choose a gauge 
for the background field A. In the following, the singular 


Bp{x) = Ap{x) + ap{x), 
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gauge is chosen to the non-perturbative instanton field 
configurations as 


Af^aix) = —r]a^iu{.X - z)^(j){x - z), 
9s 


with 


(j){x — z) = 


P 


{x - zY[{x - zY + p^y 
and the corresponding field strength tensor is 


9s _ 


{x-z)^{x-z)p 1 


^ Vaiyp 


{x — zY 


P - -5 


{{x-zY + p^y 


( 11 ) 

( 12 ) 


(13) 


with z and p denote respectively the center and size of 
the instanton, called collective coordinates together with 
the color orientation, and r]afj,i> is the ’t Hooft symbol 
which should be replaced with the anti-’t Hooft one rjapi, 
for an anti-instanton field. For the sake of simplicity, in 
the practice the most used is the spike size distribution 
n(p) = nS{p — p), where n is the overall instanton den¬ 
sity and p is the average instanton size. The fact that 
the strong coupling constant ps emerging in the denom¬ 
inator of the rhs of (ED reveals the nonperturbative na¬ 
ture of these classical configurations. In fact, instantons 
play quite important role in QCD sum rule. Early QCD 
sum rules neglecting instanton-induced continuum contri¬ 
butions didn’t obtain reliable results in many cases, but 
they are then solved by including such instanton-induced 
effects [3011^ . 

Before starting with the contraction between the quan¬ 
tum fields, we note that the time-development of the in¬ 
stanton vacuum produces the pre-exponential factor for 
the distribution of the instantons |52l[TSll53) . and is 

understood as taking ensemble average over the collective 
coordinates besides taking the usual vacuum expectation 
value due to the separation m 


^ [ dpn{p) [ dS(12|T{O,(x)Ot(0)}(”‘)|12), (14) 
ij 


where the super index ’(int)’ indicates the corresponding 
quantity containing only the interference part between the 
quantum and classical ones. Using the spike distribution, 
m becomes 


7T(“‘)(a;) = 2n J d^z{n\T{Op{x)Ol{Q)Y^^^'>\n), (15) 


with = —iDy^. Keeping only terms proportional to 
F, one has[5S] 


J d^xe^‘^-^Vyi{x,0) =^xiv-z) 

{y ~ Yp^pa-{Y9(T ^ f ^ \ 

- ^9s - ^4 - 0 fj,^{z) H-I 


2 


where the first term in rhs of the above equation is the 
pure-gluon propagator in the usual Feynman gauge, and 
the second and third ones are the leading contribution 
of the instanton field to the gluon propagator. For short 
distance region, we assume that the contribution from a 
single instanton is dominant over multi-instantons |57) . At 
the leading loop level, the gluon propagator ED becomes 
the pure-gluon one. 

In calculation, we expand the current Op into terms 
which are the products of quantum gluon fields and their 
derivatives with coefficients being composed of the instan¬ 
ton fields 


1 q 

Op{pc) — Oi{x) (1^) 

i=l 

where the operators Oi in terms of instanton and quan¬ 
tum gluon fields are listed in Appendix [C] Eq. (ITKll can be 
rewritten as 

x^/dS/ d^xF^-^{n\T{OYx)OjY))}W) 
= f;7Tf*^(g2) + .-- (19) 

where the • • • denotes the contributions from the products 
of operators being proportional to , and the expressions 
of in terms of Oi are shown in Appendix [P] 

The corresponding twelve kinds of Feynman diagrams as 
shown in the FIG. [U where the contributions from the 
first three diagrams are of the order of Og, and the contri¬ 
butions of the remainders are superficially of the order of 
al, and those from the diagrams (4),(5) and (6), in fact, 
are vanishing because of violating the conservation of the 
color-charge, namely 

yj0nt)(^2) ^ 0, fori = 4,5,6 (20) 


where the factor 2 comes from the mutually equal contri¬ 
butions of both instanton and anti-instanton. Next impor¬ 
tant step is to specify the form of the gluon propagator 
which in the background field Feynman gauge can be read 
from the part of S')!!] quadratic in a [54ll55] 

Vf,{x,y) = {0\T{al{x)al{y)}\0) 



Now, we are in the position to evaluate the contribu¬ 
tions of the remainder diagrams in FIG.[T] Using the stan¬ 
dard technique to regularizing the ultraviolet divergence 
in the modified minimal subtraction scheme, the result for 
the interference part of the correlation function is 

7T™*((5^) = cottgfiTr + aln[ci + C2{Qp)~‘^ 

+ [c3((5p)^+C 4-bC5((5p)“^] In^l , (21) 
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Fig. 1. Feynman diagrams for the interference contribution 

^{int)(g2) 

up to order a^, where spiral lines, dotted lines and 
the lines with circles denote gluons, instantons and the instan- 
ton field strength tenser respectively, and cross stands for the 
position of instantons. 


where we have ignored terms being proportional to the 
positive powers of which vanish after Borel transforma¬ 
tion, and the dimensionless coefficients Ci are numerically 
determined to be 

Co = -118.23, Cl = -3700.59as, C2 = -2394.470^, 

C3 = 11561.90as, C4 = 1850.300*, cg = 1197.24a*. (22) 

through a tedious calculation. It should be noted that 
there is no infrared divergence as expected by the instan- 
ton size being fixed in the liquid instanton vacuum model. 
Comparing Eq. m with our previous result [371138] . they 
differ not only in some coefficients but also in the loga¬ 
rithms structures due to the fact that the newly improved 
calculation is free of infrared divergence while our old one 
were not, and it would need a corresponding cutoff to reg¬ 
ularize the integral in the infrared limit. 

Putting everything above together, our final correla¬ 
tion function for the pseudoscalar glueball current is of 


the form 


yjQCD(g2^ 


-2^TT'^ny'^Kl{y) 
agfi [coTT -I- Cl -I- C2y~‘^ 

{csy^ + C4 + cgy"^) In ^ 



Oq + 


ai 




(23) 


with y = Qp. 

Before going on, let us compare the roles of the vari¬ 
ous parts of the contributions in correlation function. The 
imaginary part of the correlation function (1231) can be 
worked out to be 


ilmiT‘^®°(s) = J2{py/s)Y2{p\/s) 

TT 

+aln [czp^s - C4 -I- c^[p^s)~^] 




ao + 2ai In ■ 


-I- ( 3 In — - TT^ ) 02 


(24) 


where the pure classical contribution (the first term on 
rhs of (l24l) l is most dominant, and the contribution of the 
interference terms (the second term on rhs) is of the sec¬ 
ond place, the pure perturbative contribution simply plays 
the role of the third place, as shown in FIG. [2] where the 
imaginary part of the correlation function is multiplied 
with a weight function exp(—(s — s)/4r) as required by 
the Gaussian sum rules, and in accordance with the spirit 
of the semiclassical expansion. All these three contribu¬ 
tions are positively definite as expected. We note that the 
contribution from the topological charge screening, from 
(I5il) . is displayed in FIG. [5] as well, and its role is almost 
insignificant. Moreover, it is easy to see from FIG. [2] that 
the imaginary part of the correlation function is already 
positive from s = 0.5GeV^ to s = lOGeV^ without includ¬ 
ing the contribution of the topological charge screening, 
and the so-called positivity problem is no longer there. 
Therefore, the interference contribution to the correlation 
function is significant important not only in its magnitnde 
but also in restoring the positivity to the spectral function. 

We note here that the so-called condensate contribu¬ 
tion to the correlation function of the pseudoscalar glue- 
ball current is proven to be very small in comparing with 
the one of (1^5)) , as shown in Appendix [E] where the com¬ 
parison between the real and imaginary parts of the cor¬ 
relation function and condensate contribution to it are 
made, and shown in FIG. [5]and[51 For the reasons given 
above, the contributions from the topological charge screen¬ 
ing effect and the usual condensates are omitted in our 
sum rule analysis. 
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gluon world, is a linear combination of three operators 
J|, jg and Op after renormalization. Therefore, one as¬ 
sumes that there may be some isosinglet quark-antiquark 
pseudoscalar states mixed with the pseudoscalar glueball 
ground state G. 

Now we construct the spectral function for the corre¬ 
lation function of the pseudoscalar glueball current. The 
usual lowest one resonance plus a continuum model is used 
to saturate the phenomenological spectral function: 

ilm7T™™(s) =p“^°(s)-6»(s-so)-Im7T‘5^°(s), (33) 

TT TT 


Fig. 2. The contributions to the imaginary part of the corre¬ 
lation function from the pure instanton (cross line), interfer¬ 
ence (dashed-dotted line), pure perturbative (dotted line) and 
topological charge screening (dashed line) and the total con¬ 
tribution without the topological charge screening one (solid 
line) versus s. 


3 Spectral function 

In the isosinglet channel there are five gauge-invariant 
composite operators with the quantum numbers of 0 ^ 
which are bilinear in the fundamental quark, antiquark 
and gluon fields, namely the pseudoscalar quark densities, 
the divergences of the axial quark currents and the gluon 
anomaly: 

t8,0 


Jr = *g75(A«’72)<?, 
rJli = nrt.i5{r°/2)q], 


(25) 

(26) 
(27) 

where A® is the flavor Gell-Mann matrix, and A° = 
with I being the 3x3 flavor unit matrix. Only three of 


Op — ex. 




where sq is the QCD-hadron duality threshold, 9{s — sq) 
the step function and p^^^(s) the spectral function for 
the lowest pseudoscalar glueball state. In the usual zero- 
width approximation, the spectral function for a single 
resonance is assumed to be 




(34) 


where m is the mass of the lowest glueball, and F is the 
coupling constant of the current to the glueball defined as 


The threshold behavior for 


(0|Op(0)|G)=K (35) 

(s) is known to be 

^HAD(^) for s 0 (36) 

from the low-energy theorem in the world of no light quark 
flavors [55] or the one in the world with three light flavors 
and mu,d “C ms[59]. In fact, the threshold behavior (l36)) is 
only proven to be valid near by the chiral limit; it may not 
be extrapolated far away. Therefore, instead of considering 


these operators are independent due to the two renormalizati^fe-® coupling F as a constant [7], we choose a model for F 
invariant axial Ward identities 


as 


8 F 

Op 0^5 


7 

= 

1 


+ md + ^ma)Jl 


F = 


/ AqS, 

1 fiT, 


for s < m. 
for s > m: 


(37) 


-|--(4m„ -I- 4m£; -|- 2ms) Jl 


(28) 


B P 


72 

3 


( to „ + md- 2ms) Ji 


where the Aq and / are some constants determined late in 
numerical simulation. 

To go beyond the zero-width approximation, in fac¬ 
ing the near-actual situation, the Breit-Wigner form for a 
single resonance is assumed for p^^'^(s) 


2 

+ -^{mu +rnd- 2ms) J° + 


73 

47271 


Op. (29) 


HAD 


(s) = 


F^mF 


(s — -I- F^/4)2 + m'^F'^ ’ 


(38) 


Further, under renormalization and the flavor space rota¬ 
tion, we have 

{mJ^’°)(^r) = (30) 

= dpjp^, {dfiJpO(r) = ZdpJp^, (31) 

(07^ = 4 + (32) 

where the quantities with subscript (r) are the renormal¬ 
ized ones. As a consequence, the gluon anomaly opera¬ 
tor Op even through renormalization-invariant in the pure 


where F is the width of the lowest glubeball. 

Further, the one isolated lowest resonance assumption 
is questioned from the admixture with quarkonium states, 
and it is known from the experimental data that there 
are hve 0 pseudoscalar resonances till and around the 
mass scale of 1.405 GeV (namely 77(548), 77 ( 958 ), 77(1295), 
77 ( 1405 ) and 77 ( 1475 )). The form of the spectral function 
for hve resonances is taken to be 


= E 


is 


FfmiF, 


mlFFr^Y + mlFr 


(39) 
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where rrii and Fi being the mass and width of the i-th 
resonance, respectively. For the sake of simplicity, all cou¬ 
pling constants Fi for s < are hxed with the same Aq 
as shown in (EZl). 

It is noticed that there are other pseudoscalar reso¬ 
nances ?7(1760) and 77(2225) which are omitted from the 
summary table of PDG. These two resonances are ex¬ 
cluded in our consideration. The reasons may be listed in 
order: First, although 77(1760) and ?7(2225) may certainly 
be coupled to the pseudoscalar glueball current via the 
gluon anomaly, such coupling, however, contains a factor 
of the running coupling as commonly seen in QCD, and 
becomes weaker with increases, as demonstrated in an 
effective QCD low energy theory m- Second, the mixing 
between the considered pseudoscalar glueball and 77(1760) 
and 77 ( 2225 ) is believed to be very small because the loca¬ 
tions of 77 ( 1760 ) and 77 ( 2225 ) are far away from the scale 
of the lowest pseudoscalar glueball. Third, the continuum 
threshold sq, determined in our sum rule approach, is only 
in an effective sense due to the accuracy level of the present 
calculation. 


for the phenomenological contributions to the sum rules, 
and 


^QCD 


(so,s,r) 


g2''°is,r)-g 


CONT 

k 


(so,s,t). 


(45) 


for the theoretical contributions, where °^'^(s 0 ) s, t) is 
the contribution of continuum being defined as 


VdTTT Jso 4t 

and g^^^{s,T) is defined as 


ImQ^°(s) 

TT 

(46) 


^QCD 




2t ^ r (s + fQ^)fe77QCD(g + 7Q2) 
V47rT . 

(s-7Q2)fei7QCD(-_ -^2)- 

iQ^ 


(47) 


4 Finite width Gaussian sum rules 

In this section, we construct the appropriate sum rules 
of 0 ^ pseudoscalar glueball current which has the form 
due to dispersion relation 




(40) 


where ImTT (s) could be simulated by the phenomenolog¬ 
ical one Im7T^^^'^(s) whithin an assumed model of the 
spectral function in a spirit of sum rule approach. Using 
the Borel transformation m 


B = lim {QY 

7V-4-00 I (iV - 1)! 


dQ^ 


N 


(41) 


to both sides of OT, a family of Gaussian sum rules can 
be formed to be m 


xHAD 


(so,s,t) = 5^°°(so,s,t) 


■^=exp[-^]7T(0)4-1, (42) 
v47rr 4r 


Substituting the correlation function (^5)) of 0 + pseu¬ 
doscalar glueball into (SZl), one can derive the Gaussian 
sum rules oi k = —1, 0 and -|-1 


^QCD(~^ t) = J • 16Tr^np'^J2 (p'/s) ^2 (pY) s 
1 


TITT Ckg Cq 


nat 


'J Fkt 


exp 


s 

4t 


, —.-(ci-b C37)exp 

V47rr f 

-b C2p^\/^e“®^/®'^£>_i(-s/v^) 
- (c4/p^ - 05 ( 7 - 1)/P^)^ exp 


4r 


s 

4t 


+ n- — D^2{-al\Yr) [ao - (27 - 2)ai 

v47rT 

-b 0.5(67^ - 127 - 7 r 2 )a 2 ] 2 re"""/®^, ( 48 ) 


where sq is the continuum threshold which hadronic physics 
is (locally) dual to QGD above it. 

7J(0) = (8^)2 (43) 

+ md 


comes from the low-energy theorem for QCD with three 
light flavors [ 53 , and 


^HAD 

yfe 


(so,s,t) 



dss^ exp[— 


4t 


pHAD(g) 


(44) 


r) = X • J 2 (pv^) >2 (pv^) 

-b na\ J— : (c2p^2Te“® /®^D_2(-s/V^) 
v47rT t 

— C3\/^e“'* D-li—s j\f2T) 

-b (c4/p^ - C57 /p^) exp ^ 

+ }- — D_3(-s/v^)[2ao -b (6 - 47)01 

v47rT 

-b ( 67 ^ — I 87 — TT^ -b 6 )a 2 ]( 2 r)^^^e“'’ ^*'^(49) 
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^QCD(.^ r) = J • J2 {p\fs) Y2 {p\fs) 

+ r2c2p2(2^)3/2g-^V8r^_3(_^/y^) 

y A'kt ^ 

- C 32 Te“® /®^£>_2 (-s/v^) 

H-^^=-D_ 4 (—s/\/^)[6ao + (22 — 127)01 

v47rT 

+ (187^ - 667 - Stt^ + 36)a2]4r2e"®"/®^, (50) 


where 


J = 


\/47rr , 


ds exp 


(s-s)2 


4r 


(51) 


and parabolic cylinder function D-d-iisy/r) defined as 


D^d-i{z) = 




d\ 


(dz)« 


d > 0. 

(52) 


5 Numerical simulation 


The expressions for the three-loop running coupling 
constant as{Q^) with three massless flavors {Nf = 3) at 
renormalization scale p |62] 


as{p^) 




1 

CW 


^4/ 


Po 


(53) 


are used, where 44(4^)/tt is the two-loop running cou¬ 
pling constant with {Nf = 0) 


44 jp'^) _ 1 Pi InL 

TT PoL Po {PqLY 


(54) 


and the values of the average instanton size and the overall 
instanton density are adopted from the instanton liquid 
modeipS] 


n = lfm “4 = 0 . 0016 GeV^ 

-p = ifm ~ 1.667GeV^ (57) 

O 

The resonance parameters in Ea. dMl) could be esti¬ 
mated by matching both sides of sum rules Eq. (021) op¬ 
timally in the fiducial domain. In doing so, the parame¬ 
ter s and the threshold sq should be determined priority. 
Eirstly, it is obvious that sq must be greater than the mass 
square of the highest lying isolated resonance considered, 
namely 


So > TOmax ( 58 ) 

in our multi-resonance assumption (or just the resonance 
mass itself in the case of a single resonance assumption), 
and should guarantee that there is a sum rule window 
for Gaussian sum rules. Secondly, the peak positions of 
s, t) versus s curves should not change too much 
with moderate variation of sq. Here we do not mention 
about the values of r in this condition, because the peak 
positions of these curves is not affected by appropriate val¬ 
ues of r. It is found that the behavior of these curves can 
satisfy above requirements as shown in FIG. [3] if sq lies 
in the interval of (4GeV^, 5GeV^). It is also remarkable 
that the peak positions has already indicated the approx¬ 
imate mass of the hadron considered. Thus, we would ex¬ 
pect that the mass of the pseudoscalar glueball should be 
nearby the value 4§ ~ 1.449 GeV. In another way, if one 
uses the curves of I/^^^(so, s, r) versus r with fixed s and 
So to obtain the physical parameters through 0iT|. then 
s should be set approximately to be Speak of the curves 
of t/^^^(so, s, t) versus r, so as to highlight the underly¬ 
ing hadron state in consideration and suppress the con¬ 
tributions from other states. Besides, it needs a sum rule 


and 


L = In 
1 


, Li = Im L — InL — 1, 


Po = 


Pi = 


1 

42 


11--iV, 


102 - ^Nf 


43 


2857 5033 325 

- Nf H- Nf 

2 18 ^ 54 ^ 


(55) 


with the color number N^ = 3 and the QCD renormal¬ 
ization invariant scale A = 120MeV. We take p^ = 4^ 
after calculating Borel transforms based on the renormal¬ 
ization group improvement for Gaussian sum rules [63]. 
The subtraction constant 7T(0) has been fixed as |7] 


7T(0) ~ -0.022GeV^ 



(56) 


Fig. 3. The curves of s, r) versus s with different sq 

at r = 1 GeV‘‘. 


window which the hadron physical properties should be 
stable in this region. For the upper limit r^ax of the sum 
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rule window, the resonance contribution should be great 
than the continuum one 

^^fe^°('SO,S,Ti„ax) > Tmax) (59) 


according to the standard requirement due to the fact that 
in the energy region above Tmax the perturbative contri¬ 
bution is dominant. At Tmin, which lies in the low-energy 
region, we require that the single instanton contribution 
should be relatively large so that 




> 50% 


(60) 


In the same time, to require that the multi-instanton cor¬ 
rections remain negligible, we simply adopt a rough esti¬ 
mate 

r„ta > 2p-‘. (61) 

According to the above requirements, we find that in the 
domain 

T G (0.5,4.5)GeV'^, (62) 


our sum rules work very well. Finally, in order to measure 
the compatibility between both sides of the sum rules dUl) 
in our numerical simulation, we divide the sum rule win¬ 
dow [Tmin, Tmax] hito N — 100 Segments of equal width, 
[n,Ti+i], with To = Tinin and tn = Tinax, and introduce a 
variation 5 (called the matching measure) which is defined 
as 


_ 1 ^ [L{n)-R{n)Y 
\L{r^)R{n)\ 


(63) 


where L(ri) and R{Ti) are Ihs and rhs of evaluated 
at Ti. 

Let us first consider the case of single-resonance plus 
continuum model p4l) of the spectral function by exclud¬ 
ing the interference contribution i7™*((5^) from 
(the case A), in order to recover the results before. In 
this case, the imaginary part of {Q"^), however, be¬ 

comes negative for s below 3.9GeV^, and the full interac¬ 
tion plays a role of repulsive potential in the pseudoscalar 
channel. This is the reason why the authors in Ref. [39] 
cannot find the signal for the pseudoscalar glueball. When 
we chose to work in the positively definite region, say 
s G (4,10)GeV^, The mass of the 0 glueball can be 
worked out by using the family of Gaussian sum rules 
dHj). The fitting parameters are listed in the first three 
lines of Tab. [Tjand the corresponding matching curves for 
fc = — 1, 0 and -1-1 are displayed in FIG. S] respectively. 
The optical values of mass, coupling constant and sq of 
0 pseudoscalar glueball are 


m = 2.010 ± 0.299GeV, / = 0.584 ± 0.043GeV, 

So = 5.21 ± 0.43GeV^ (64) 

where the errors are estimated from the uncertainties of 
the spread between the individual sum rules, and by vary¬ 
ing value of A in the region of (the same for hereafter) 

(65) 


The mass values of 0 + pseudoscalar glueball in (IMll are 
reasonably consist with the one obtained in Ref. H by 
adding the topological charge screening effect and per¬ 
forming the so-called Gaussian-tail distribution for instan¬ 
ton size. However, after performing the Gaussian-tail dis¬ 
tribution for instanton size, the mass of the 0++ glueball is 
lower to be around 1.25GeV [7] in contradiction with the 
lattice simulation and phenomenology [65] . The 

mass scales in (IMl) locate at the strong repulsive potential 
region of energy (below 3.9GeV^) where the bound state 
of glueball cannot form when working in the spike distri¬ 
bution which is motivated from the liquid instanton model 
of QGD vacuum in the large W limit. In fact, the funda¬ 
mental spectral positivity bound can be traced back to 
the definition of the correlation function (jSj). The spectral 
function should be positive even before taking the average 
with any specific instanton-size distribution. It is difficult 
for us to understand that there is no artificial in chang¬ 
ing the positivity behavior of the spectral function just by 
performing the average with the Gaussian-tail distribution 
for the instanton size. 

From now on, the full correlation function (Q^) 

including the interference contribution 77™*(Q^), (l23l) . is 
used in our analysis of the Gaussian sum rules (H^ . In the 
case of single-resonance plus continuum models, specified 
respectively by (IMI) and (l38ll . for the spectral function, 
the optimal parameters governing the sum rules with zero 
(the case B) and finite (the case G) widths are listed from 
the fourth to ninth line of Tab. [T] and the corresponding 
curves for the Ihs and rhs of (HD) with /c = — 1, 0 and 
-1-1 are displayed in FIG. [5] and [5] respectively. From Tab. 
m the optical values of the pseudoscalar glueball mass, 
width, coupling and the duality threshold with the best 
matching are: 

m = 1.644 ± 0.194GeV, / = 1.412 ± 0.129GeV, 

So =4.77±0.74GeV^ (66) 

for one zero-width resonance model, and 

m = 1.407 ± 0.162 GeV, R = 0.053 ± 0.018 GeV 

/ = 1.687 ± 0.145 GeV, sq = 4.63 ± 0.62 GeV^ (67) 

for one finite-width resonance model. It is shown in FIG. 
[S] that the topological charge screening effect has little 
impact on Gaussian sum rules indeed. 

In the numerical simulation for the case of the five 
finite-width resonances plus continuum model (1391) for the 
spectral function (the case D), we just choose the data in 
PDG as the fitting parameters for masses and width of 
the resonances ?7(548), ?7(985), ?7(1295) and 7^(1475), and 
the result of the single resonance model (the case G) as 
the fitting parameters for 77(1405); while the couplings of 
the five resonances to the current are chosen to be ap¬ 
proximately the same as that for 77(1405) determined in 
case C because the pseudoscalar qurkonia can be directly 
coupled to the gluon anomaly, and as a consequence, all 
five resonances should be coupled to the current with the 
strengths of almost the same magnitude of degree; finally, 
the optimal parameters are determined by adjusting the 


A = 120 - 200MeV. 
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chosen parameters so that the matching measure <5 for 
both sides of the Gaussian sum rules (1421) is minimal. The 
optimal parameters governing the sum rules are listed in 
the remaining lines of Tab. [1] The corresponding curves 
for the Ihs and rhs of (ga) with k = —1, 0 and +1 are 
displayed in the FIG. [71 Taking the average, the optical 
values of the widths of the five lowest 0 ^ resonances in 
the world of QCD with three massless quarks, and the 
corresponding optical fit parameters are predicted to be 

w^( 548 ) = 0.548 ± 0.022 GeV, /^( 548 ) = 1.133 ± 0.167 GeV, 



TlGeV^ 


71,(548) = 1-3 X 10"® ± 3.9 X 10“® GeV (68) 


mn{ 985 ) = 0.958 ± 0.051 GeV, /^(ggs) = 1.200 ± 0.233 GeV, 
71,(985) = 1-9 X 10-3 ± 5.7 X 10-3 GeV (69) 

"*,,( 1295 ) = 1-295 ± 0.075 GeV, /^(i 295 ) = 1-202 ± 0.112 GeV, 
71,(1295) = 0.055 ± 0.018 GeV (70) 

^r,(i 405 ) = 1.405 ±0.081 GeV,/^(i 405 ) = 1.313 ±0.105 GeV, 



71,(1405) = 0.051 ± 0.017 GeV (71) 

w^(i 475 ) = 1.475 ± 0.092 GeV,/^(i 475 ) = 1.023 ± 0.097 GeV, 
71,(1475) = 0.085 ± 0.028 GeV (72) 

with 

So = 4.78 ± 0.64 GeV^ (73) 

The FIG. [6] and [7] show the satisfactory compatibility be¬ 
tween both sides of the sum rules over the whole fiducial 
region. These results are in good accordance with the ex¬ 
perimental discovered resonance |14] 

771^(1405) = 1409.8 ± 2.5MeV, 71^(i405) = 51.1 ± 3.4MeV 

(74) 



Fig. 4. The Ihs (dashed line) and rhs (solid line) of the sum 
rules dill) with k = —1, 0,1 versus r in the case where the cor¬ 
relation function contains only pure perturbative 

contribution and pure instanton one, and a single zero-width 
resonance plus continuum model is adopted for the spectral 
function. 


6 Discussion and Conclusion 

The instanton-gluon interference and its role in finite- 
width Gaussian sum rules for 0 pseudoscalar glueball 
are analyzed in this paper. Our main results can be sum¬ 
marized as follows: 

First, the contribution to the correlation function aris¬ 
ing from the interference between the classical instanton 
fields and the quantum gluon ones is reexamined and de¬ 
rived in the framework of the semi-classical expansion of 
the instanton liquid vacuum model of QGD. The resul¬ 
tant expression is gauge invariant, and free of the infrared 
divergence, and differs from our previous one not only in 
some coefficients but also in the logarithms structures [571 
155] . Its magnitude is just between the larger contribu¬ 
tion from pure classical instanton configurations and the 
smaller one from the pure quantum fields, and plays a 
great role in sum rule analysis in accordance with the spirit 


of semi-classical expansion. The imaginary part of the cor¬ 
relation function including this interference contribution 
turns to be positive without including of the topological 
charge screening effect which is proven to be smaller than 
the perturbative contribution in the fiducial sum rule win¬ 
dow, and negligible in comparison with the interference 
effect. The so-called problem of positivity violations in 
the imaginary part of the correlation function, stressed 
by H. Forkel[7], disappears. Moreover, it is excluded in 
the correlation function the traditional condensate con¬ 
tribution to avoid the double counting[7] because conden¬ 
sates can be reproduced by the instanton distributions [53 
mUTlIlH]; Another cause to do so is that the usual con¬ 
densate contribution is proven to be unusually weak, and 
cannot fully reflect the nonperturbative nature of the low- 
lying gluonia[551[71 [571[55] : In our opinion, the condensate 
contribution can be considered as a small fraction of the 
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Table 1. The optimal fitting values of the mass m, width F, coupling constant /, continuum threshold so , s, and matching 
measure 5 for the possible 0 ^ resonances in the sum rule window [rmin, Tmax] for the best matching between Ihs and rhs 
of the sum rules (mj with k = —1,0,1 are listed, where in case A, the correlation function contains only pure 

perturbative contribution and pure instanton one, and a single zero-width resonance plus continuum model is adopted for the 
spectral function, while all the contributions arising from pure instanton, pure perturbative and interference between both are 
included in the correlation function for cases B, C and D, in which a single zero-width resonance plus continuum model of the 
spectral function is adopted for case B, and a single Hnite-width resonances plus continuum model for case C, and the five 
hnite-width resonances plus continuum model for case D, respectively. 


case 

k 

resonances 

^(GeV) 

m(GeV) 

r(GeV) 

/(GeV) 

so(GeV^) 

[rmin, 'rmax](GeV"‘) 

S/IQ-^ 


-1 


2.050 

1.950 ± 0.050 

0 

0.561 

5.26 ±0.14 

[0.8, 4.0] 

22.60 ± 1.46 

A 

0 


2.100 

2.000 ± 0.035 

0 

0.599 

5.18 ±0.12 

[0.5, 6.0] 

1.99 ±0.13 


1 


2.110 

2.080 ± 0.040 

0 

0.592 

5.20 ±0.10 

[0.8, 6.0] 

7.56 ± 0.45 


-1 


1.405 

1.682 ±0.023 

0 

1.457 

5.24 ±0.15 

[1.0, 3.0] 

49.72 ± 2.98 

B 

0 


1.490 

1.620 ±0.027 

0 

1.387 

4.61 ±0.09 

[0.8, 2.5] 

3.67 ±0.18 


1 


1.380 

1.631 ±0.031 

0 

1.392 

4.45 ±0.11 

[1.4, 3.2] 

8.97 ± 0.43 


-1 


1.405 

1.405 ± 0.024 

0.05 

1.630 

5.25 ±0.14 

[0.8, 4.0] 

2.44 ±0.12 

C 

0 


1.350 

1.400 ± 0.025 

0.08 

1.671 

4.45 ±0.12 

[0.5, 4.3] 

4.47 ±0.21 


1 


1.380 

1.416 ±0.024 

0.03 

1.760 

4.19 ±0.11 

[0.5, 4.5] 

2.19 ±0.11 



77(548) 


0.548 ± 0.008 

1.3 X 10”'’ 

1.100 






77(958) 


0.958 ±0.014 

1.9 X 10"® 

1.100 





-1 

77 ( 1295 ) 

1.405 

1.295 ±0.020 

0.055 

1.200 

5.25 ±0.12 

[0.5, 4.5] 

4.88 ± 0.23 



77(1405) 


1.405 ±0.021 

0.051 

1.330 






77 ( 1475 ) 


1.475 ±0.023 

0.085 

1.010 






77 ( 548 ) 


0.548 ± 0.009 

1.3 X 10"° 

1.200 






77 ( 958 ) 


0.958 ±0.016 

1.9 X 10"® 

1.300 





0 

77 ( 1295 ) 

1.500 

1.295 ±0.021 

0.055 

1.195 

4.79 ±0.11 

[0.5, 4.3] 

3.38 ±0.17 

D 


77(1405) 


1.405 ± 0.022 

0.051 

1.300 





77 ( 1475 ) 


1.475 ± 0.025 

0.085 

1.011 






77 ( 548 ) 


0.548 ±0.010 

1.3 X lO"'’ 

1.100 






77 ( 958 ) 


0.958 ±0.017 

1.9 X 10"® 

1.200 





1 

77 ( 1295 ) 

1.405 

1.295 ±0.022 

0.055 

1.210 

4.30 ±0.10 

[1.0, 4.0] 

3.56 ±0.17 



77(1405) 


1.405 ± 0.023 

0.051 

1.310 






77 ( 1475 ) 


1.475 ± 0.024 

0.085 

1.050 





corresponding instanton one, so it is naturally taken into 
account already. 

Second, the properties of the lowest lying 0 ^ pseu¬ 
doscalar glueball are systematically investigated in a fam¬ 
ily of Gaussian sum rules in five different cases. In case 
A, the correlation function [Q"^) contains only pure 

perturbative contribution and pure instanton one, and 
a single zero-width resonance plus continuum model is 
adopted for the spectral function, and of course, the old re¬ 
sults are recovered (even excluding the topological charge 
screening contribution), and some pathology is explored. 
To go beyond the above constraint, all the contributions 
arising from pure instanton, pure perturbative and inter¬ 
ference between both are included in the correlation func¬ 
tion for cases B, C and D, in which a single zero-width 
resonance plus continuum model of the spectral function 
is adopted for case B, and a single finite-width resonances 
plus continuum model for case C, and the five finite-width 
resonances plus continuum model for case D, respectively. 
The optimal fitting values of the mass to, width F, cou¬ 
pling constant /, continuum threshold sq for the possible 
0 ^ resonances are obtained, and quite consistent with 


each other. The main difference between this work and 
our previous one [571155] is that for the spectral function of 
considered resonances, instead of the zero-width approx¬ 
imation of one gluonic resonance plus another low-lying 
quark-antiquartk ones, the finite-width Breit-Wigner form 
with a correct threshold behavior for the lowest five res¬ 
onances with the same quantum numbers is used in this 
work in order to compare with the phenomenology. The re¬ 
sultant Gaussian sum rules with A: = — 1,0, -1-1 are carried 
out with a few of the QCD standard inputting parameters, 
and really in accordance with the experimental data. 

As a discussion, let us now identify where the low¬ 
est lying 0 ^ pseudoscalar glueball is. The result of the 
single-resonance plus continuum models B and C, namely 
Eqs. (1551) and (1571) . imply that the meson 77(1405) may be 
the most fevered candidate for the lowest lying 0 ^ pseu¬ 
doscalar glueball because the difference between the two 
models is just the width of the resonances, and the latter 
is of course believed to be more in accordance with the re¬ 
ality. This conclusion can further be justified by the result 
of the five-resonances plus continuum model, namely Eqs. 
(ED, (ED and (ED- Note that the first two resonances 
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T [GeV] 




T[GeV^ 


Fig. 5. The Ihs without the topological charge screening contri¬ 
bution (dashed line), rhs with the topological charge screening 
contribution (dot line) and rhs (solid line) of the sum rules (1 42 II 
with k = —1,0,1 versus r in the case where the interference 
contribution is included in the correlation function 
and a single zero-width resonance plus continuum model is 
adopted for the spectral function. 


?7(548) and 77 ( 986 ) are far away from the mass scale of 
77 ( 1405 ), and usually considered as the superposition of the 
fundamental flavor-singlet and octet pseudoscalar mesons 
composed of quark-antiquark pair to have a dominant role 
in responsible for the axial anomaly [571[551[55] . In order to 
explore the structures of the remainder three resonance 
77 ( 1295 ), 77 ( 1405 ) and 77 ( 1475 ), we would like use the rj-rj'-G 
mixing formalism based on the anomalous Ward identity 
for transition matrix elements |491[T^ to relate the physical 
states 77 , 77 ' and G to the fundamental flavor-singlet and 
octet quark-antiquark mesons and the lowest pure gluon 





Fig. 6. The Ihs (dashed line) and rhs (solid line) of the sum 
rules (O with k = —1,0,1 versus r in the case where the 
correlation function (Q^) contains the pure instanton, 

interference, and pure perturbative contributions, and a single 
finite-width resonance plus continuum model is adopted for the 
spectral function. 

state through a rotation 

/(0|Op|?7i295)\ /(0|Op|778)\ 

(0|6p|77i475) = (OlOplm) . (75) 

V(0|Opbi405)/ V(0|Op|G)/ 

where the U is the mixing matrix mM 

( cos (fip — sin (fp \ 

cos(f>GSm(pp cos (^p sin (/)G sint/fc (76) 
— cos ipG sin (pp — cos (pp sin (pa cos pc I 
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-QCD 
' HAD 


Ml = 0.548GeV fi = l.lOOGeV 

ri = l,3xlO-®GeV 

M 2 = 0,958GeV /2 = l.lOOGeV 

r2 = 1.9xlO-^GeV 

M 3 = l,295GeV /a = 1.200GeV 

Fa =0,055GeV 

Mi = l,405GeV fi = 1.330GeV 
r4 =0.051GeV 

Ms = 1.475GeV fs = l.OlOGeV 
r5 =0.085GeV 


T [GeV^] 




after normalization, respectively. As a relatively rough es¬ 
timation, from Eq. (f75|l the values of the couplings of Op 
to the states 771 , rjs and pseudoscalar glueball G are ob¬ 
tained 


(0|6p|?7i) = -0.117GeV^ 

(0|6p|?78) = 0.568GeV^ 

(0|dp|G) = 0.813GeV^ (79) 

by reversing Eq. Ga). This shows that the coupling of 
the 0“+ glueball current to the pure glueball state G is 
dominant, and the signs of the couplings ( 0 |Op| 77 i) and 
(OlOpIrys) are similar to those predicted by the scalar glueball- 
meson coupling theorems [701158]. 

Furthermore, 77(1405) and 77(1475) could originate from 
a single pole (see Amsler and Masoni’s review for eta(1405) 
in PDG). To check whether the single pseudoscalar meson 
assumption may be consistent with our sum rule approach, 
or the dependence of the results on the model selected for 
the spectral function, we add an analysis of four resonance 
model for the spectral function. The result is shown in Ap¬ 
pendix F. From FIG. ITOland Tab. 2 where it is easy to see 
that the four resonance model for the spectral function 
is inconsistent with our sum rule analysis because, firstly 
the resultant mass scale of the fourth resonance is approx¬ 
imately 1.41GeV in average which locates outsides of the 
sum rule window; secondly the match degree S of the four 
resonance model is obviously worse than the 6 of the five 
resonance model. 

In summary, our result suggests that 77(1405) is a good 
candidate for the lowest 0 “+ pseudoscalar glueball with 
some mixture with the nearby excited isovector and isoscalar 
qq mesons. This is a first theoretical support for the phe¬ 
nomenological estimation from the sum rule approach. 


Acknowledgements 


Fig. 7. The Ihs (dashed line) and rhs (solid line) of the sum 
rules (|4^ with k = —1,0, 1 versus r in the case where the 
correlation function (Q^) contains the pure instanton, in¬ 

terference, and pure perturbative contributions, and five finite- 
width resonance plus continuum model is adopted for the spec¬ 
tral function. 

where (pp « 40° and (j)c ~ 22° is 77 — 77 ' [JS] mixing angle 
and the mixing angle of the pseudoscalar glueball with 77 '. 
Then, one has 

/ 0.766044 -0.642788 0 \ 

U = 0.595982 0.286965 0.374607 (77) 

y-0.595982 -0.286965 0.927184 J 

To be quantitative, the corresponding normalized cou¬ 
plings F to the three resonances 77 , 77 ' and G with masses 
1.295, 1.475, 1.405 GeV can be read from Tab. [1] to be 

0.51GeV^ 0.65GeV^ 0.56GeV^ (78) 


This work is supported by BEPC National Laboratory 
Project R&D and BES Gollaboration Research Founda¬ 
tion. 

A Pure instanton and perturbative 
contribution 

In instanton liquid model, the pure instanton contri¬ 
bution of 0 ^ pseudoscalar glueball with spike instanton 
distribution in Euclidean space-time is known as[7i r7Tir72l 

[75] 

= -2^n^np‘^Q^Kl{Qp). (80) 

where iLl is the McDonald functions and n{p) is the size 
distribution of instantons, and the appearance of the over¬ 
all minus sign in rhs of (l80l) is due to the anti-hermitian 
property of the current, i.e. = —Op in Euclidean space- 
time. 
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The second part of cni) has already been calculated up 
to three-loop level in the MS dimensional regularization 
scheme [TTlFT^rTi] 


C The operators Oi{x) 

The operators Oi in terms of instanton and quantum 
gluon fields are 


77(qu)(Q2) ^ 


Q' r ^ 1 Q" 

tto -I- Oi in — 


^ 1 2 Q" 

-I- 02 in — 


(81) 


where /x is the renormalization scale, and the coefficients 
with the inclusion of the correct threshold effect are 


1 -h 20.75 -h 305.95 


oo — —2 

oi = 2 ^ 

02 = -10.1250 


f^s\ 

9 /as\ 

— 

- -f 72.531 — 

V TT / 

4 V TT /_ 


(82) 


Oi{x) = F^y^a[A{x)]Fpa^a[A{x)] 

02{x) = ^Fpy^a[A{x)]{dpaaa[A{x)]) 

0^(^X^ — ‘^9sfabcFp,i/^a[-^{.x')^Api)(^X^Q,(jc{,X^ 

04 {x) = 4:{dp,a,,a{x))idpacra{x)) 

0^{x^ — ^QsfahcApl){x^{dp^(Xi/a{xy)(lac{x) 

0^{x^ — ^9 Q fabcfadeAp^l){x^Ap(i{^X^(Xi/c{x^(X(je.{x^ 

0'j{x^ — ^9sfabcFp^jj^a\A{x^\(Xpl){x^(Xiyc{x^ 

^ 8 (^) — ^9sfabcO>pb{x^(Xi/c{x^{dp(X(ja{xy) 

Oq{x^ — ^9 s f(^bcfadeAp(i{x^Clpl){x^(Xi/c{x^^<yeix^ 

^io(^) — 9 s f(^bcfade0^iib{x^(lyc{x^Qjp(i{F)(l(yQ{x^ 

where Fp,ij^a[A{x)\ is the instanton field strength associ¬ 
ated with the instanton field A. 


for QCD with three massless quark flavors up to three- D The interference contributions n'r^\q‘) 
loop level. 

The expressions of the interference contributions [q^) 
in terms of Oi are 


B Topological charge screening 


The topological charge screening effect can be under¬ 
stood by tracing back to the anomaly of the axial-vector 
current = T'iV'iTsT/it/’i with 'ijji being quark field of fla¬ 
vor i 


d^Jl = 2NfQ{x) (83) 


where the topology charge current Q{x) relates to the 
pseudoscalar glueball current as Q{x) = Op{x)/[Sir). Eq. 
(|83ll indicates that instantons generate quark-antiquark 
pairs, which then may form a light meson to mediate the 
long-range multi-instanton interaction m- Kikuchi and 
Wudka [7S] suggest that such a light meson, which can 
couple to the instanton in the way that the instanton gen¬ 
erates all light quark flavors with the identical probability, 
be the meson 770 , and constructs an effective Lagrangian 
(see also [7^177] ! which gives rise to the topological charge 
screening contribution 77 *°p to the correlation function of 
the pseudoscalar current as [7] 


7T‘°p(Q2) 


V 




+ 


T 


-I- 


(84) 


where rriri and are the masses of the mesons 77 and 
77 ', and the two constants and F^i are evaluated to be 
Wirn^sincj) and Wirn^ cos cj), and ^ and (j) ~ 22 ° are the 
coupling strength of rjo to an instanton and the rj — rj' 
mixing angle, respectively. 


where 


yj(ci+qu)(^ 2 ) ^ y^^|O 2 (a:)O 2 ( 0 )|f 2 ) 
Tjf+‘^")(g2) = f(f2|O2(x)O3(0)|f2) 
Tjf+‘^"^(< 72 ) = f(f2|O3(:r)O3(0)|f2) 
yjM+qu)(^ 2 ) ^ f{n\Oiix)05m^} 
yj(cl+qu)(^ 2 ) ^ y^^|O 4 (^)Og( 0 )|| 2 ) 
Tjf+‘^"^(< 72 ) = f{n\Oiix)07m^} 
7 j(cl+qu)(^ 2 ) ^ y^^|O^(^)( 35 ( 0 )|| 2 ) 
Tjf+‘^")(g 2 ) = f{[2\Orix)Ormn) 
Tjf+‘^"^(( 72 ) = f(f2|O5(:r)O7(0)|f2) 
= f{n\Oeix)07m^} 
n[f+^'^\q^) = f{n\05ix)Oemn) 
n[f+^'^\q^) = f{n\Oeix)Oem^} 

T — gjicp^i/pfjCpf y' pf(j' f d z f d X. 


( 86 ) 

(87) 


E Comparison between the condensate 
contributions and the instanton-induced ones 

The contributions arising from condensates up to the 
eighth dimensions are known to be as follows [7] 

^co„d(Q2) ^ + lal{asG^)\n% 

7T 

- S>a^{gG^)— + —^al{asG'^)'^-^, (88) 
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where 


(a.G^) = O.OSGeV^ 

{gG^) = 0.27GeV^(asG^), (89) 

The imaginary part of condensates (1551) has the form 


F The results for the four finite-width 
resonance model from the Gaussian sum rules 

The optimal parameters of the four finite-width reso¬ 
nance model for the spectral function are listed in Tab [5J 
and the corresponding plots are shown in FIG. ITTll 


—Im7T^°“'^(s) = ——a^(asG^) — 8a^{gG^)S{s) 

TT TT 

-I— —al{asG'^)^S'{s). (90) 

The comparison between the imaginary part of the cor¬ 
relation function, (IMll . and condensate contribution to it 
are shown in FIG. [H while the comparison between the 
various real parts of (1^ (which altogether are related 
with the imaginary part by the dispersion relation (PO])) 
and condensate contribution to it are shown in FIG. M 



Fig. 8. The contributions to the imaginary part of the corre¬ 
lation function from the condensate (dashed-dotted line) and 
the total contribution (solid line) versus s. 



Fig. 9. The contributions to the correlation function aris¬ 
ing from the pure instanton (cross line), interference (dashed- 
dotted line), pure perturbative (dotted line) and topological 
charge screening (dashed line) and the total contribution with¬ 
out the topological charge screening one (solid line) versus Q^. 





Fig. 10. The Ihs (dashed line) and rhs (solid line) of the sum 
rules (liSl) with k = —1,0,1 versus r in the case where the 
correlation function (Q^) contains the pure instanton, 

interference, and pure perturbative contributions, and the four 
finite-width resonance plus continuum model is adopted for the 
spectral function. 
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